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The combination of a vortex line in a one-dimensional optical lattice with fermions bound to 
the vortex core makes up an ultracold superstring. We give a detailed derivation of the way to 
make this supersymmetric string in the laboratory. In particular, we discuss the presence of a 
fermionic bound state in the vortex core and the tuning of the laser beams needed to achieve 
supersymmetry. Moreover, we discuss experimental consequences of supersymmetry and identify 
the precise supersymmetry in the problem. Finally, we make the mathematical connection with 
string theory. 
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I. INTRODUCTION 

Ultracold quantum gases provide a very exciting 
branch of physics. Besides the interesting physics that 
the gases offer by themselves, it has also been possible 
in the last few years to model with quantum gases sys- 
tems from other branches of physics, and by doing so to 
provide answers to long-standing questions. The latter 
is mainly due to the amazing accuracy by which their 
properties can be tuned and manipulated. This involves 
the trapping potential, the dimensionality, the interac- 
tion between the atoms, and the statistics. By using a 
three-dimensional optical lattice the supcrfluid-Mott in- 
sulator transition in the Bosc-Hubbard model has been 
observed [l| . Bosonic atoms confined in one-dimensional 
tubes by means of a two-dimensional optical lattice where 
shown to realize the Lieb-Liniger gas 0, ■ The uni- 
tarity regime of strong interactions was reached by us- 
ingFeshbach resonances to control the scattering length 
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To this shortlist of examples from condensed-matter 
theory, also examples from high-energy physics can be 
added. In a spinor Bose-Einstein condensate with fer- 
romagnetic interactions skyrmion physics has been stud- 
ied 0, H3j whereas an antifcrromagnctic spinor Bose- 
Einstein condensate allows for monopole or hedgehog so- 
lutions [ll], [lj ■ There is also a proposal for studying 
charge fractionalization in one dimensio n [la, and for 
creating (static) non-abelian gauge fields [lj, [l5j . In re- 
cent work [l|j we have added another proposal to model a 
system from high-energy physics. By combining a vortex 
line in a one-dimensional optical lattice with a fermionic 
gas bound to the vortex core, it is possible to tune the 
laser parameters such that a nonrelativistic supersym- 
metric string is created. This we called the ultracold su- 
perstring. This proposal combines three topics that have 
attracted a lot of attention in the area of ultracold atomic 
gases. These topics are vortic es Il7l \MJ 03 . I20L l2lj . Bose- 
Fermi mixtures |22l I2I I2I I2I I2tll2ll2if7and optical lat- 
tices [j], I23 ■ Apart from its potential to experimentally 
probe certain aspects of superstring theory, this proposal 
is also very interesting because it brings supersymmetry 



within experimental reach. 

Supersymmetry is a very special symmetry, that re- 
lates fermions and bosons with each other. It plays an 
important role in string theory, where supersymmetry 
is an essential ingredient to make a consistent theory 
without the so-called tachyon, i.e., a particle that has 
a negative mass squared. In the physics of the mini- 
mally extended standard model, supersymmetry is used 
to remove quadratic divergences. This results in a super 
partner for each of the known particles of the standard 
model. However, supersymmetry is manifestly broken 
in our world and none of these superpartners have been 
observed. A third field where supersymmetry plays a 
role is in modeling disorder and chaos [3(|. Here super- 
symmetry is introduced artificially to properly perform 
the average over disorder. Finally, supersymmetry plays 
an important role in the field of supersymmetric quan- 
tum mechanics, where the formal structure of a super- 
symmetric theory is applied to derive exact results. In 
particular this means that a supersymmetry generator 
Q is defined, such that the hamiltonian can be written 
as H = {Q,Q^}, which is one of the basic relations in 
the relativistic superalgebra. It is important for our pur- 
poses to note, that this relation is no longer enforced 
by the superalgebra in the nonrelativistic limit. Careful 
analysis |3lL l32 | shows that in this limit the hamiltonian 
is replaced by the number operator, i.e., J\f = {Q,Q^}- 
It may sometimes be possible to write a nonrelativistic 
hamiltonian as the anticommutator of the supersymme- 
try generators, but this does not correspond to the non- 
relativistic limit of a relativistic theory 

In our proposal, a physical effect of supersymmetry is 
that the stability of the superstring against spiraling out 
of the gas is exceptionally large, because the damping 
of the center-of-mass motion is reduced by a destruc- 
tive interference between processes that create two ad- 
ditional bosonic excitations of the superstring and pro- 
cesses that produce an additional particle-hole pair of 
fermions. Moreover, this system allows for the study of 
a quantum phase transition that spontaneously breaks 
supersymmetry as we will show. 

Another very interesting aspect of the ultracold super- 



string is the close relation with string-bit models |33| . 
These are models that discretize the string in the spatial 
direction, cither to pcrturbativcly solve string theory, or, 
more radically, to reveal a more fundamental theory that 
underlies superstring theory. String-bit models describe 
the transverse degrees of freedom of the string in a very 
similar fashion as in our theory of the ultracold super- 
string. 

In this article we investigate in detail the physics of 
ultracold superstrings, expanding on our previous work 
[lq . The article is organized as follows. In Sec. II we 
give the detailed derivation of the conditions for the ul- 
tracold superstring to be created. In particular, we pay 
attention to the presence of the fcrmionic bound state in 
the vortex core and the tuning of the lasers to reach su- 
persymmetry. In Sec. Ill we investigate the experimental 
consequences of the supersymmetry. Sec. IV contains a 
detailed description of the supersymmetry by studying 
the supcralgcbra. In Sec. V we make connection with 
string theory Finally, we end with our conclusions in 
Sec. VI. 




FIG. 1: (Color online) Artist's impression of the setup. The 
disks represent the bosonic condensate density and the blue 
balls in the vortex core represent the fermionic density. The 
black line is a guide to the eye to see the wiggling of the vortex 
line that corresponds to a Kelvin mode. 



persymmetry is possible. The temperature is taken to be 
well below the Bosc-Einstcin condensation temperature, 
so that thermal fluctuations are strongly suppressed. We 
only consider the zero-temperature limit, because super- 
symmetry is formally broken for nonzero temperatures. 



II. ULTRACOLD SUPERSTRINGS 



Our proposal makes use of the fact that a vortex line 
through a Bose-Einstein condensate in a one-dimensional 
optical lattice can behave according to the laws of quan- 
tum mechanics y4|. Such an optical lattice consists of 
two identical counter-propagating laser beams and pro- 
vides a periodic potential for atoms. When applied along 
the symmetry axis of a cigar-shaped condensate, which 
we call the z axis from now on, the optical lattice divides 
the condensate into weakly-coupled pancake-shaped con- 
densates. In the case of a red-detuned lattice, the gaus- 
sian profile of the laser beam provides also the desired 
trapping in the radial direction. Rotation of the Bosc- 
Einstcin condensate along the z axis creates a vortex line 
that passes through each pancake. Quantum fluctuations 
of the vortex position are greatly enhanced in this con- 
figuration because of the small number of atoms Nb in 
each pancake, which can be as low as Nb = 10, but 
is typically around Nb = 1000. An added advantage 
of the stacked-pancake configuration, as opposed to the 
bulk situation, is that the dispersion of the vortex oscilla- 
tions is particle like. This ultimately allows for supersym- 
metry with the fermionic atoms in the mixture. In the 
one-dimensional optical lattice the vortex line becomes 
a chain of so-called pancake vortices. This produces a 
setup which is pictured schematically in Fig ^ There 
is a critical external rotation frequency f2 c above which 
a vortex in the center of the condensate is stable. For 
O < tt c the vortex is unstable, but because of its Euler 
dynamics, it takes a relatively long time before it spirals 
out of the gas [l^, yj, |3g| . We analyze in detail the case 
of fi = 0, i.e., the situation in which the condensate is no 
longer rotated externally after a vortex is created. How- 
ever, the physics is very similar for all £1 < fi c , where su- 



A. Atomic species 



A convenient choice for the boson-fermion mixture is 
87 Rb and 40 K, since such Bose-Fermi mixtures have re- 
cently been realized in the laboratory J2J|, |2J, J25L l26j . 
and because the resonance lines in these two atomic 
species lie very nearby. The mostly used \f, m/) hy- 
perfine spin states are |9/2,±9/2) and |9/2,±7/2) for 
40 K, and |2,±2), |1,±1), and |1,0) for 87 Rb J37J. They 
all have a negative interspecies scattering lenght asF, 
which is not desirable for our purposes as we show be- 
low. It could be possible to use other spin states, which 
have a positive interspecies scattering length. An other 
possiblity is to tune the scattering length, using one of 
the various broad Fcshbach resonances that can make 
the interaction repulsive while keeping the probability to 
create molecules negligible |37j . 

In principle it is also possible to use other mixtures. 
Another Bose-Fermi mixture that has been realized in the 
laboratory consists of 23 Na and 6 Li atoms [22|,|27]]. This 
mixture is less convenient because the resonance lines 
are widely separated, so that the two species feel very 
different optical potentials and it is hard to trap both 
with a single laser. In addition, 6 Li is relatively hard to 
trap in an optical lattice because of its small mass. For 
these reasons, the 23 Na- 6 Li mixture can only be used in a 
very restricted parameter regime, as we will show lateron 
in Fig. \5\ For the same reasons, the mixture 87 Rb - 
6 Li J28I does not work well either. The mixture 7 Li- 6 Li 
[38|,|33 cannot be used at all, because the resonance lines 
of the species are the same, so it is impossible to tune the 
physical properties of the mixture. 
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FIG. 2: (a) (Color online) Schematic picture of the setup. 
Here r is the radial distance in the xy plane. The pink and 
blue blobs represent the bosonic and fermionic densities, re- 
spectively. Moreover, A is the wavelength of the laser. The 
blue and red lines indicate the strength of the optical poten- 
tial, respectively, for the bosons and fermions as a function of 
the z coordinate, (b) Schematic fine structure level scheme 
of the bosonic and fermionic atomic species. Because we con- 
sider only sufficiently large detunings the hyperfine level struc- 
ture is not resolved. 



B. Optical lattice 

Because the excited states of the bosonic and fermionic 
atoms have different transition frequencies, the optical 
lattice produces for the two species a periodic potential 
with the same lattice spacing, but with a different height, 
as schematically shown in Fig. [21 This is very crucial, be- 
cause it allows to tune the optical lattice for the bosonic 
and fermionic atoms scpcratcly, by careful adjustcmcnt 
of the wavelength and the Rabi frequency, i.e., the inten- 
sity of the laser. This is required to be able to tune the 
system to become supersymmetric lateron. 

For the 86 Rb- 40 K mixture the Rabi frequencies are 
in a good approximation the same. For other mixtures 
the Rabi frequencies are different and we then take the 
bosonic Rabi frequency as a reference. We take into ac- 
count the fine-structure level scheme of the atoms, but, 
assuming that we are sufficiently far from resonance, we 
neglect the hyperfine structure. As a result, the optical 
potential is given by 



V b ,f(z) = V b ,f cos 2 (2irz/\), 



where the well depths obey 



V b ,f = - 




(1) 



(2) 



uj = 2itc/\ is the laser frequency, and lod 1 and u>d 2 arc 
the frequencies of the D\ and D^ resonance lines. Here we 
neglected spontaneous emission of photons. This effect 
slightly modifies the trapping potential, but gives a finite 
lifetime to the atoms. Using the rotating-wave approx- 
imation and neglecting the fine structure, the effective 



rate of photon absorption can for red-detuned laser light 
be estimed as 
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where Tg j? is the lincwidth of the bosonic or fermionic 
excited state, respectively. For blue-detuned laser light, 
the atoms are trapped in the regions of low laser intensity 
and spontaneous emission is strongly reduced. 

The optical potential should be sufficiently deep to 
have a bound state for the bosonic and fermionic atoms. 
To make sure that that is the case we impose the condi- 
tion 
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where we have used the recoil-energy 
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which is the energy associated with the absorption of a 
photon. On the other hand, the optical lattice should not 
be so stron g to drive the system in the Mott-insulator 
state U I29L |4CJ | . In one dimension with many atoms per 
site, this requires an exceptionally deep lattice, which 
only occurs if the laser frequency is very close to the 
resonance frequency of the atomic species. Since we stay 
away from resonance, this situation does not occur in our 
calculations. 

The wavefunctions in the z direction arc assumed to be 
the groundstate wavefunctions of the harmonic oscillator 
associated with the optical lattice and thus given by 
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For the tunneling amplitude, we use the expression [4 
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which becomes exact for a deep lattice. Therefore, the 
atomic dispersions along the z axis are given by 



es.F(fc) = 2J b ,f[1 - cos(fcA/2)]. 



(9) 



Lateron we need for the fermions the relation between 
the average number of particle per site and the chemical 
potential fiF- From the above dispersion we derive at 
zero temperature that 
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(10) 



where we neglect also interaction effects. 



C. Kelvons 

The wavefunction in the (axial) z direction is fully 
specified by the optical lattice and all the dynamics 
thus takes place in the radial direction, i.e., in the xy- 
plane. Since the vortex-fluctuations form the lowest-lying 
modes, we restrict the dynamics to the vortex motion. 
We follow the derivations in earlier work J2J, |42j , where 
a specific ansatz for the wavefunction was used, to achieve 
this. In this work the condensate density was described 
by a gaussian wavefunction with size -Rtf and the vortex 
core was approximated by a step function. Furthermore, 
it was assumed that the vortex is close to the center. The 
motion of the vortices results in kelvons, i.e., quantized 
oscillations of the vortex, described by the creation and 
annihilation operators 



-Rtf 



(x + iy), ir 
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(x-iy), (11) 



which obey [b, b'] = 1. Without the optical lattice Kelvin 
waves have already been observed [43L lil ] . The kelvons 
have the dispersion 
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bosonic harmonic length in the radial direction, and u> 
is the associated frequency. Using another ansatz for the 
condensate wavefunction can slightly change the constant 
of proportionality in the definition of the kelvon operators 
and in the details of the dispersion, but the dispersion al- 
ways stays tight-binding like. 

For the calculation of the bound state in the vortex 
core, we need to go beyond the description of the core 
by a step function. This change of the calculation could 
improve the value of Jk, but not the functional form 
of the kelvon dispersion. Since the corrections on the 
value of Jk are small, we just use the result in Eq. Ijl3|) . 
Besides the bandwidth Jk we derive from Eq. IJ12J) also 
the chemical potential for the kelvons, which gives 
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Note that the chemical potential is positive only for suf- 
ficiently small rotations, which is due to the fact that 
the vortex is in principle unstable for these values of the 
rotation and wants to spiral out of the center of the gas 
cloud. 

D. Bound states in the vortex core 



where 



j K = r[o,(e/R TF ) 4 ]j B , 



(13) 



r[0, z] is the incomplete Gamma function, -Rtf is the 
Thomas-Fermi radius in the radial direction, £ is the 



By treating the interaction between the bosonic and 
fcrmionic atoms in mean-field approximation, we have 
to solve the Gross-Pitaevskii equation for the condensate 
wavefunction \I/(r) coupled to the Schrodingcr equation 
for the fermion wavefunction ip{ T ) 
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which we investigate for the case that fl = 0. The inter- 
action paramaters are related to the scattering lengths 
according to 
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with clbb the boson-boson scattering length and <xbf 
the boson-fcrmion scattering length and m,R the reduced 
mass rnBTUp / {viib + %)■ Although it is very well pos- 
sible to solve these equations numerically, we prefer an 



analytic treatment, to gain more insight into the problem. 
To proceed we make the approximation that the conden- 
sate wavefunction is not affected by the presence of the 
fermions. This is justified, because the contribution of 
the fermions is Np/Ns smaller than the contribution of 
the bosons, where Nb,f is the average number of bosons 
and fermions at a lattice site. This ratio will be smaller 
than 10 -3 as it turns out. Taking into account the inter- 
action with the fermions leads to a slightly wider vortex 
core, which enhances the possiblity of a bound state. So 
we first solve the Gross-Pitaevskii equation for the con- 
densate density neglecting the presence of the fermions 
and then use the condensate density as an effective po- 



tcntial for the fcrmions. Since we only want to estimate 
when there is a bound state and we do not need the 
details of this bound state, we make the following ap- 
proximations. First, we assume the vortex to be in the 
center such that the problem is rotationally symmetric 
and we only have to solve the radial equation. Because 
of the quantum uncertainty the vortex position in prin- 
ciple fluctuates around the center of the trap, but these 
fluctuations are small. Second, we assume the envelope 
condensate wavefunction to be Thomas-Fermi like, i.e., 



If we take for the healing length £ the usual expression 
in the center of the trap, i.e., 
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we obtain the relation 
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Third, we describe the vortex core by r 2 /(r 2 + 2£ 2 ) |45j . 
such that the total bosonic density is given by 
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By expressing the energy in terms of /ibi^b /mp we can 
write the Schrodingcr equation for the fcrmions as 
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where e = m F E/ fi B mB and the dimensionless parameter 
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determines whether or not there is a bound state in the 
core of the vortex. 

If we assume that £ <C £ and £ <C if , we can neglect 
the harmonic confinement and the Thomas-Fermi profile 
of the Bose-Einstein condensate. The effective potential 
for the fermionic atoms is then given by Tr 2 /(r 2 + 2£ 2 ). 
This potential has a bound state for each value of T, 
because for large distances from the core, it behaves as 
r(l — 2£ 2 /r 2 ). However, the size of the wavefunction 
describing the bound state becomes extremely large for 
small values of I\ Hence it is necessary to take into ac- 
count the exact form of the potential to make a quanti- 
tative estimate of the existence of the bound state. The 
potential is determined by the values of the radial bosonic 
and fermionic harmonic length £ and £ F . Since i F deter- 
mines the potential outside the condensate, it determines 
whether or not the fermions can tunnel out of the core 
to this region. For the existence of the bound state we 
can neglect this contribution, which is always justified, 
because it enhances the possibility of having a bound 
state. 



The radial bosonic harmonic length i is fixed by the 
normalization of the condenstate wavefunction 



■n TF (r) = N B . (25) 
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Neglecting the presence of the core we find the usual 
expression for the Thomas-Fermi profile 
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we see that taking into account the core implies that we 
have to solve the equation 

t_ Am B N B U BB 

£ 4 TTh 2 

(27) 

where the last two terms come form the presence of the 
core. Since the core is small in this approximation, this 
results in a radial harmonic length that is only slightly 
modified. The requirement that the wavefunction should 
vanish well within the condensate can then be quantified 
to yield the expression 

a 2 e* 



a 2 + 2£ 2 
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(28) 



FIG. 3: (Color online) Setup with the additional laser to ma- 
nipulate the radial trapping potential. 



where <r is the radial size of the fermionic wavefunction. 
In this way we obtain that for typical densities there is a 
bound state for T > 1.5 which means that ubb/cibf > 2. 
In contrast to the radial bosonic length t, the fermionic 
radial harmonic length t r F is not fixed. When the optical 
lattice is red-detuned, the lattice can be used to trap the 
atoms also in the radial direction. As a consequence, the 
total confining potential for the fermions is a multiple of 
the confining potential of the bosons, i.e., 
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This gives the relation 
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from which we derive 
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However, if the lattice is blue-detuned or if the radial 
trapping is tuned independently, this relation is not true. 
The radial trapping can be tuned by introducing a second 
running laser in the same direction as the optical lattice, 
as shown in Fig [31 The new laser beam has a constant 
intensity along the z axis, and does not influence the 
one-dimensional potential wells, but it does change the 
radial confinement. In principle this second laser also 
introduces interference terms, but they are much faster 
than the atoms can follow for the frequencies of interest 
to us. Therefore, the intensities of the two lasers can sim- 
ply be added. In particular, as we show lateron, adjusting 
the radial trapping potentials is needed to get supersym- 
metric interaction terms. The condition imposed by this 
requirement is 
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which gives the following expression for the fermionic ra- 
dial harmonic length 
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FIG. 4: Effective potential for the fermions. Lengths are 
measured in units of £ and energies in units of HBmB/rriF- 
(a) r = .5, Pp 3> £: No bound state, since the potential is 
too small, (b) F = 2.6, t r F S> I: Bound state in the core, 
but possibility to tunnel outside, (c) F = 2.6, £ r F ~ £: Bound 
state in the core, no tunneling possible. 



In this last case, the harmonic radial potential for the 
fermions is very small. In principle this allows the 
fermionic atoms to tunnel out of the vortex core, to the 
region where the condensate density vanishes. However, 
the tunneling is suppressed by increasing the parameter 
r. A WKB estimate gives that for V > 5 the lifetime 
of the fermions in the core is larger than a second. This 
means that asB /o-bf > 8. Further increasing this ratio 
increases this lifetime dramatically. Since adjusting the 
radial trapping potentials is only needed close to the cen- 
ter of the trap, it is also a possibility to use a second laser 
with a much smaller waist, such that higher-order con- 
tributions from the potential prevent the fermions from 
tunneling out of the core. For various situations, the 
effective potential for the fermions is shown in Fig0J 



E. Interactions 



In our supcrstring realization there are also boson- 
boson and boson-fermion interactions. The kelvons in- 
teract repulsively among each other when il < fl c . For 
fl = the kelvon-kelvon interaction coefficient is given 



by |4C 
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In addition, a repulsive interaction between the kelvons 
and the fermionic atoms is generated by the fact that 
physically the presence of a kelvon means that the vor- 
tex core is shifted off center, together with the fermions 
trapped in it. Because of the radial confinement experi- 
enced by the trapped fermions, this increases the energy 
of the vortex. When the vortex core is shifted from to 
r, the fcrmion hamiltonian is extended by a term 
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where e'e is the number operator for the fermions in the 
core. Defining Cb,f as the spring constants associated 
with the radial confinement of the bosonic and fermionic 
atoms, respectively, and using the definition of the kelvon 
operators, this translates into 



Z1\b 



(35) 



So the kelvon-fermion interaction coefficient is found to 
be 
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FIG. 5: Tuning of the lattice laser to obtain supersymmetry. 
Plotted is the Rabi frequency for the bosonic atoms versus 
wavelength for 87 Rb- 40 K for 10000 (solid line), 1000 (dashed 
line) and 500 (dotted line) bosonic atoms per site. Note that 
for the blue-detuned part, i.e., A < 760 nm for the 87 Rb- 40 K 
mixture extra radial trapping is needed, either magnetically, 
or by using an extra running laser as discussed in the text and 
shown in Fig. |3] In Fig. |5]we display how to tune the running 
laser to obtain also supersymmetric interactions. In the inset 
we plotted the Rabi frequency that is required for the 23 Na- 
Li mixture to obtain supersymmetry, again for 10000 (solid 
line), 1000 (dashed line) and 500 (dotted line) bosonic atoms 
per site. Note that this can only be obtained in a very limited 
range of wavelength's. 



F. Supersymmetry 

To obtain a supersymmetric situation we have three 
requirements. In the first place the hopping amplitudes 
have to be the same 



Jf = Jk = t. 



(37) 



This can be done by adjusting the laser parameters A 
and VLb, as shown in Fig[5] The freedom in choosing the 
wavelength of the laser can be used to minimize the atom 
loss. In Fig. El we plot the atom loss as a function of the 
wavelength of the laser. 

Secondly, the chemical potentials have to be the same 



\x F = \x K = £i. 



(38) 



This can be achieved by adjusting the fermion filling frac- 
tion Np, as shown in Fig[7| Using the result from Eq. 
(|1U[) and using the requirements for supersymmetry we 
obtain 



Np = — arcsin 

7T 



--"TT7 



0. 






\ 4j B r 



0, 
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: — arcsm 

7T 



V 
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z b\ m(v B /E B y 



/4p 



0. 



The ratio 1/1% is undetermined by supersymmetry con- 
straints. In order for the Thomas-Fermi approximation 
to apply in the radial direction, versus the gaussian wave- 
function in z direction, this ratio needs to be sufficiently 
small. In the figure a ratio of 1/5 is chosen. 

Finally, the interaction terms have to be the same. 
This implies 



Vkk = Vkf = U. 



(39) 



Setting these coefficients equal to each other gives a con- 
dition on the radial trapping given by 



Cf_ 
C b 



Rtf 



r[o, {i/Rtf) 



(40) 



The radial trapping can be tuned by introducing a sec- 
ond running laser, as explained before. For the second 
laser, we can again independently choose both the wave- 
length and the Rabi frequency as shown in Fig. [H] This 
can again be used to minimize the atom loss due to the 
red-detuned laser, but it turns out that atom loss is al- 
ways quite small anyway for reasonable system parame- 
ters. Only for very small detunings, the lifetime is less 
than a second. 



0.1 




- 


0.08 




- 


0.06 




- 


0.04 




- 


0.02 


j ^ 







700 800 900 

Wavelength (nm) 



S 5 











§0.08 

|0.04 

\ 1 
X. 4 




,^j , , ~"~~ 


/^- 


DO 


600 800 1000 1200/ 
Wavelength (nm) / 


- 


\^ 



600 700 800 900 1000 

Wavelength running laser beam (nm) 



FIG. 6: Effective rate of photon absorption as a function of 
the wavelength of the optical lattice for 1000 bosonic atoms 
per site. The solid line is for the bosonic atoms, whereas 
the dash 
mixture. 
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FIG. 7: Tuning of the average number of fermions per lat- 
tice site to obtain supersymmetry for 10000 (solid line) , 1000 
(dashed line) and 500 (dotted line) bosonic atoms per lattice 
site. The result depends on the ratio of the bosonic harmonic 
lengths in the axial and radial directions £/£%■ This ratio 
should be sufficiently small to be radially in the Thomas- 
Fermi limit. For this plot a ratio of 1/5 is chosen. 



FIG. 8: Tuning of the additional laser to obtain supersymmet- 
ric interactions for 1000 bosonic atoms per lattice site. The 
Rabi frequency is plotted versus the wavelength of the running 
laser for different wavelengths of the lattice laser beam: 1000 
nm for the left curve and 600 nm for the right curve. In the 
inset the effective rate of photon absorption originating from 
the additional laser is plotted, again the left curves are for a 
lattice laser of 1000 nm and the right curves are for a lattice 
laser of 1000 nm. The solid lines indicate the rate of photon 
absorption for the bosonic atoms, whereas the dashed lines 
are the rates of photon absorption for the fermionic atoms. 



Here bi is the annihilation operator of a kelvon at site 
i, Ci is the annihilation operator of a fermion at site i, 
(ij) means that the summation runs over neighbouring 
sites, and fj! = /x — 2i. We used the convention for the 
Fourier transformation fk — (l/-\/N s ) J2 n elkZn fn, where 
N s is the number of lattice sites. We define a = A/2 
as the lattice spacing and L = N s a as the length of the 
system. Assuming that N s 3> 1, such that L > a, wc 
can perform a continuum approximation to obtain for the 
hamiltonian 



n = 



+ / dzct(z) ("^t~^) c(z) 



(42) 



U 



+- / dz [tf{z)b\z)b(z)b{z) + 2b\z)b{z)c\z)c(z)}, 



G. Hamiltonian 

Combining everything, our superstring is described by 
the supersymmctric hamiltonian 



n = -*5>t&j + 4< 



(41) 



(U) 



-Y\-^{b\b l + cU) 



jb]b\bibi + Ub\b lC \ Cl 



where we introduced the effective mass m* = h 2 /2a 2 t. 
This continuum hamiltonian turns out to be exactly solv- 
able [43, E3 by a straightforward generalization of the 
Bethe-Ansatz solution of the one-dimensional Bose gas 
|49l l50j . However, the exact solutions spontaneously 
break supersymmetry and do not give much insight in 
the role of supersymmetry in the problem. 
Using that the lagrangian is given by 



u*4»' 



d 



c}ih-^-c. 



H[b\b;< 



(43) 



the action in the continuum limit is obtained as 



H. Nonlocal interaction 



,5' 



dt \ dz { b* \ih— 



if- 



2m* dz 2 
h 2 d 2 



+ fi)b (44) 



H c 



1 dt ' 2m* dz 2 



which now explicitly shows the supersymmetry of the 
problem, because it remains invariant when b and c are 
rotated into each other. If we neglect the interaction 
terms, which are rather small anyway, the fermions fill 
a Fermi sea and the low-energy excitations are particle- 
hole excitations around the Fermi surface. Therefore, 
the low-energy part of the theory is properly described 
by linearizing the fcrmionic dispersion around the Fermi 
level. To preserve supersymmetry we do the same for the 
bosons and obtain at the quadratic level the action 

S= J2 I dt I ' dkh* a {k,t)(ihd t - hv F {ak - k F ))b a (k,t) 

<T=±-' J 

+c* (fc, t)(ihdt — hvp(ak — k F ))c a (k, t) >, 

where a indicates whether the particles are right 
movers or left movers and v F = hk F /m* is the 
Fermi velocity. We used that fi = Y?k 2 F j2m* . 
We identify the Fermi velocity with the veloc- 
ity of light c and perform the transformation 
c a {z,t) = (l/\/I) £fe ^^(A; + ak F ,t). We intro- 
duce the Dirac spinor i[>(z,t) = (c+(z,t),c_(z,t)) and 
V>(z,i) = V^(-z,£)7°, with 7° = a y . The other Dirac ma- 
trices are 7 1 = ia x and 7 s = a z . The two bosonic fields 
can be captured in a single Klein-Gordon field X(z, t) = 
(l/VL)J2k etkz VWk(b+{k + k F ,t) + b-{k-k F ,t)), 
such that d^d^X = (pile 2 - d 2 )X = 0. This enables us 
to rewrite the linearized action as 



d 2 x(d fi X*d ll X 



+ itv4>^d^}, (45) 



which is the action for the transverse modes of a free rel- 
ativistic N = 1 superstring in 3 + 1 dimensions J51J . In 
modern language, the Lorentz invariance of this action 
appears here as an emergent phenomenon at long wave- 
lcnghts, because the underlying theory is not Lorentz 
invariant. This is very similar with the way in which 
Lorentz invariance appears in string- bit models |33j . A 
second property of this action is, that the fermionic 
part has classically chiral symmetry, but quantum- 
mechaniclly suffers from a chiral anomaly. Whereas in 
string theory this is an unwanted feature, in our case it 
has a physical origin, because it comes about from the 
fact that the underlying microscopic theory does not con- 
serve the chiral current itpj^^ip, and only conserves the 
current ip^/^ip associated with the conservation of the to- 
tal number of fermions. 



The presence of a kelvon implies that neighbouring vor- 
tex cores are slightly shifted with respect to each other. 
This effect decreases the fermionic hopping amplitude 
and results in a interaction term that couples fermions 
on neighbouring sites of the form 



ri! = t s £ J A(k')cos(k)c\c k b\ l b kl . 



(46) 



7v/v' 



Since this term breaks supersymmetry, we want to in- 
vestigate the system parameters for which it can be ne- 
glected, i.e., for which A(k') <C 1. To do so we consider 
a kelvon with a certain wavenumber k. The relative dis- 
tance between neighboring cores kan then be estimated 
to be 

Sr = ^kx (47) 



From Eq. 112011 we know that for small distances the vor- 
tex core can be modeled as a harmonic potential with 
width £. Hence, the fermionic wavefunctions are gaus- 
sians with the same width. So we have to compute 



A(k) 



J d 2 re 



2« 2 e 



J d 2 re 



(48) 



= 1-e 



,Vr 
4p 






fc A -ZTrprji 



l6N b l 4 ' 

where we used the relation from Eq. 126fl . From this 
same relation we see that (-Rtf/^) 4 scales with the num- 
ber of bosonic atoms Nb, such that A(k) is indepen- 
dent of Nb- We can estimate B^ f /Nb£ 4 to be of order 
unity, such that the requirement for A(k) to be small 
only depends on the wavenumber k. If we identify this 
wavenumber with the Fermi momentum, i.e., 



k 2 F \ 2 
16 



«1, 



(49) 



we obtain a restriction on the fermionic filling fraction 
which can be estimated to be 

N F < 0.1. (50) 

From Fig. \7\we sec that for most of the parameter space 
this condition is fullfillcd. 



III. EXPERIMENTAL SIGNAL 

It is an important question how the supersymmetry 
can be observed. Therefore we need to distinguish be- 
tween the question whether the hamiltonian is tuned to 
be supersymmetric and whether the quantum ground 
state is supersymmetric, since it is possible that the 
ground state can spontaneously break supersymmetry. 
We are primarely interested in the situation that both 
the hamiltonian and the quantum ground state are su- 
persymmetric. 
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FIG. 9: Diagrams that contribute to the dissipation of the 
vortex line. The solid lines represent bosonic propagators, 
the dashed lines represent fermionic propagators. The left dia- 
gram is called XlnK{k, u>), the left diagram is called TlKF(k, u>) 



of the deviation from supersymmctry. We can extend 
this to higher order correlation functions. The condition 
that 



{{b\hf) = ((4 Cl ) 2 ) = (4a), 



(53) 



can be used to prove that in order to have supersymmctry 
also the condition 



A. Density measurements 

The two observables that are most easy to measure 
experimentally are the average number of fermions at a 
site Np and the average number of kelvons Njc- The 
average fermion number can be determined by usual 
absorpsion measurements. The number of kelvons can 
be obtained from the mean-square displacement (r 2 ) = 
(1/N S ) J2i( x l+yi) °f the pancake vortices, which can be 
measured by imaging along the z direction the size of the 
circle within which the vortex positions are concentrated 
Hal- Because 



b% 



N F 



<JJ 



R 2 



1 

2' 



(51) 



this can directly be translated to the number of kelvons 
at a site. It is clear that in order to have a supersym- 
metric state, the kelvon and fermion modes should have 
the same average occupation number, i.e.. 



(b}bi) = (ch 



(52) 



This allows us to devise an experimental measure for the 
proximity to the supcrsymmctric point, which can be di- 
rectly measured, namely 



N B 



p2 

-Tlrpp 



- N F - 1/2 



This quantity has an absolute minimum of zero at the su- 
persymmetric point, so that it's magnitude is a measure 



N 



B p 4 



R 



2N f 



(54) 



IT 



should hold. The quantiy (r 4 ) can again be measured 
from the distribution of the measured vortex postions. 



B. Dissipation 

Another consequence of supersymmetry that can be di- 
rectly measured is the reduced dissipation. Dissipation in 
this context results in the vortex spiraling out of the gas. 
|35l |3j| . The dominant part of the dissipation is given 
by the coupling to the kelvon modes and the fermionic 
modes. The lowest order diagrams are given in Fig. [U] 
and denoted by Hkk (k, to) for the coupling to the kelvon 
modes and IlK F {k,uj) when there is also coupling to the 
fermionic modes. The imaginary part of these diagrams 
measures the dissipation. In order to be able to know the 
dissipation away from the supersymmetric point we per- 
form the calculation for unequal dispersions tK,F{k) and 
unequal coupling constant Ukk and Ukf- We introduce 
the usual notation for the Bose-Einstein and Fcrmi-Dirac 
distribution functions 



N B [e(k)} = 



1 



e(k)/k B T _ i' 



N F [e(k)} = 



1 



e e(k)/k B T + I' 



The diagrams are then given by 



J 



IL KK (k,iuj) 
U KF (k,iuj) 



2[/ 2 



KK 



-U 



KF 



dp 

2^ 



dp f dp' [ l + N B[ei<(p)} + N B [e K {p')])N B [e K {p + p' - k)] - N B [e K {p)]N B [e K {p')] 
2n 



2k 
dp'(- 



1 



ihuj - e K (p) - ejf(p') + e K (p + p' - k) 
N B [e K {p)]+ N F [e F {p')])N F [e F {p + p' - k)} + N B [e K {p)]N F [e F {p')] 



2- 



ihuj - e K (p) - e F (p') + e F (p + p' - k) 



where the minus sign in front of the expression for 11^-^ (k,iu>) comes from the presence of the fermion loop. Note 
that due to the different combinatorial factors the diagram T1 BB comes with an extra factor of two, which is lacking 
in the case of the diagram H B f- As a result, the two diagrams do not cancel exactly at the supcrsymmctric point, as 
we claimed previously [l6j . Instead, the dissipation is reduces by a factor 2 |52j . The imaginary part of the diagrams 
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gives the following expressions 

lm[n KK (k,cu)] = 2U 2 KK J^J^{5(h^-e K (p)-e K (p') + e K (p + p'-k))x 

((l + N B [e K {p)] +N B [e K {p l )])N B [€ K {p + p l - k)] - N B [e K (p)}Nb[ck(p')}) }, 
]m[U KF (k,w)] = -U 2 KF J^J^{6(nw-eKip)-e F (p') + e F (p+p'-k))x 

((-1 - N B [e K (p)] + N F [e F (p')fjN F [e F (p+p' - k)} + N B [e K (p)]N F [e F {p')}^ . (55) 



r 



At zero temperature we have that N F [e(k)] = 8(e(k)), 
and N B [e(k)} = -6»(e(fc)) = -N F [e(k)}. Using this, we 
see that if there is supcrsymrnetry, i.e., if ei<(k) — £ F (k) 
and Ukk = Ukf, we have that 

U KK (k,u) = -2X\ KF {k,u), 

and in particular that 

lm.[U KK (k,u!) + U KF (k,w)] = -lm[U KK (k,uj)}, 

such that at zero temperature supersymmetry results in a 
dissipation rate that is only half as large as in the case of 
a ordinary vortex-line. Using these expressions, it is also 
possible to calculate the quantum dissipation at nonzero 
temperature, or when supersymmetry is broken. In par- 
ticular, when the interaction coefficients are tuned away 
from the supersymmetric point such that U BF = \2U BB 
and supersymmetry is maintained at the quadratic level, 
the dissipation exactly vanishes and the superstring is 
extremely stable in the center of the condensate. 



breaking. As a consequence of the breaking of the U{\) 
symmetry because of the Bose-Einstein condensation, the 
dispersion of the kelvon modes becomes gapless. The dis- 
persion becomes the usual Bogoliubov dispersion, which 
reads 



Hu B (k) = ^e(k)2 + 2fie(k), 



(58) 



with e(k) = h 2 k 2 /2m*. For long wavelengths this yields 
a linear behaviour. Also the fermionic modes become 
gapless. This is a result of the breaking of supersymmetry 
and this mode is called the goldstino. Because (b) = 
^J [ijXJ , the dispersion for the goldstino is given by 



hu F {k) = e(k) - (j. + \(b)f = e(fc), 



(59) 



which results in a quadratic dispersion. Clearly the 
bosonic and fermionic dispersion in Eqs. 158|) and H59|) 
are now different, which signals a nonsupersymmetric sit- 
uation. 



C. Spontaneous supersymmetry breaking 

When the hamiltonian is supersymmetric, the ground 
state still can break supersymmetry. This is the phe- 
nomenon of spontaneous supersymmetry breaking. For 
< Q c , the ultracold superstring is unstable against 
Bose-Einstein condensation of kelvons. This breaks su- 
persymmetry, because the fermionic modes cannot Bose- 
Einstein condense. Bose-Einstein condensation implies 
that the kelvon annihilation operator obtains an expec- 
tation value 



(bi) * 0. 



(56) 



From the definition of the kelvon operator we conclude 
that as a consequence 



(xf + (yf> 0. 



(57) 



This means that the vortex moves out of the center of the 
trap. Experimentally this is easy to measure. Moreover, 
by monitoring the vortex position when it moves out of 
the center of the trap, this also allows for the experi- 
mental investigation of the dynamics of supersymmetry 



IV. SUPERSYMMETRY 

In this section we review the algebra associated with 
supersymmetric field theories both in the relativistic (su- 
per Poincare algebra) and the non-relativistic limit (su- 
per Galilei algebra). We give an explicit representation 
of the super Galilei algebra in terms of the bosonic and 
fermionic operators. 



A. Super Poincare algebra 

Associated with a rclatistic field theory in D = d+ 1 di- 
mensions is the Poincare algebra, whose generators con- 
sist of a vector P M , that generates translations and an 
antisymmetric tensor J^, that generates Lorenz tran- 
sormations. The greek indices run from to d = D — 1, 
such that P° should be identified with the hamiltonian 
7i, up to a constant. The algebra is then given by 



[P»,P V \ = 
[J» u ,pp] = i( 7 fpp» 



(60) 



r fp P n) 



[J"", J'*] = i{^PJ va +if a JP u -rfoj^ -rf a J p »), 
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where rf v is the flat space Minkowski metric. When 
there is supersymmetry we can extend this to the super 
Poincarc algebra. For N — 1 supersymmetry in 1 + 1 
dimensions this involves the two-component Majorana 
spinor Q a , a = 1,2, which is the generator of supersym- 
metry transformations. The algebra is then extended to 
include also 



[P",Q a ] -- 


= 


(61) 


[P»,Qa\ -- 


= 


(62) 


[J^,Q a ] = 


= -\[r,i v UQp 


(63) 


{Qa,Qp} = 


- H*p» 


(64) 



where the 7^ arc again the Dirac matrices. We use con- 
ventions such that Q a has two real components. To make 
connection with the supersymmetry in the ultracold su- 
perstring we combine these two components in one com- 
plex supersymmetry operator 



Qi + 1Q2 



(65) 



This decomposition breaks manifest Lorentz symmetry, 
but since we are ultimately interested in the nonrelativis- 
tic limit, this is of no concern to us here. As a result we 
obtain the following algebra 



[P^Q] = 


= 




(66) 


[J^,Q} = 


= -{^Q 1 




(67) 


{Q,0 f } = 


-. p° 




(68) 


{Q,Q} = 


-- {Q\Q ] } = 


-P 1 . 


(69) 



In particular, we see that the hamiltonian P° is fixed by 
the supersymmetry generator. This is a very peculiar re- 
striction on the hamiltonian, which is only true for the 
rclativistic theory. In the nonrclativistic limit, the su- 
persymmetry decouples from the space-time translation 
symmetry as we show now. 



generators P and the hamiltonian H. Wc can now take 
the limit c — * 00 to obtain the super Galilei algebra. The 
Galilei algebra obtained in this manner has nonvanishing 
commutators 



[P,K] 
\H,K\ 



im*N 
iP. 



(74) 

(75) 



The part involving the supersymmetry becomes only 

{Q,Q t }=AA. (76) 

This defines the algebra SiQ J3JJ. As is clear, in this case 
the hamiltonian is decoupled form the supersymmetry. 
In 1 + 1 and 2 + 1 dimensions, it is sometimes possible 
to define an extended superal geb ra 1S2S, which again in- 
volves the hamiltonian |33L J54L J55| . In d = 1 this amounts 
to introducing an extra scalar supersymmetry generator 
R with the algebra: 

{Q,i?t} = -P (77) 

{R,R?} = H/2. (78) 

C. Representation 

The representation for the S\Q algebra in terms of the 
bosonic and fcrmionic operators b and c, can easily be 
found to be 



TV 



dz{b\z)b(z) + J(z)c(z)}, 



(79) 



it) 



P = / dz{b\z)d z b(z) + c*<(z)d z c{z)}, (80) 



Q 



t _ 



dz c'(z)b(z), 

dz b\z)c(z). 



(81) 

(82) 
(83) 



B. Super Galilei algebra 

The Galilei algebra can be derived as a limit of the 
Poincare algebra by performinga Inonii-Wigner contrac- 
tion [53 in the following way |3ll l32| . We write 



p u = 

P 1 = 

J 01 = 

Q = 



- (m*Afc 2 
c 

P 
cK 



■H) 



(70) 

(71) 
(72) 
(73) 



where c is the speed of light and m* denotes the mass, 
which is the same for the bosonic and fermionic degrees 
of freedom. We also defined a number operator A/", which 
counts all the particles in the system, and boost opera- 
tors K. Furthermore, wc still have the space translation 



In addition, we can thus also define 

h 



R = 



i?t = - 



dz c ! (z)d z b(z), 
dz b t (z)d z b(z). 



(84) 
(85) 



This produces 



{R, tf> = / dz {-H*)^K*) - JM^wA*)} ■ 

(86) 
which indeed is the kinetic energy part of the hamilto- 
nian. The full quadratic part of the hamiltonian can be 
expressed as 



H = i{i?, J R t }-M0,Q t }- 



(87) 
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For completeness, we mention that we can also use su- 
perspace techniques to write the hamiltonian in a mani- 
fest supersymmetric way. This involves the introduction 
of a complex superfield 

*(z,9) = e~ e ' e/2 (b(z)+9*c(z)) (88) 

m*(z,9) = e- ' '/ 2 {b*{z)+c*{z)e) 1 (89) 

where 9 is a Grassman variable such that {9,9} = 
{9*, 9*} = and {9,9*} = 0. The hamiltonian is in 
terms of the superneld given by 



n 



d9*d9 / dz 



r h 2 

1 2m* 1 v 



/*!*(*,< 



u 



* 



{z,9t}. 



(90) 



In this formulation the spontaneous breaking of super- 
symmetry is particularly elegant, because the hamilto- 
nian has the form of a standard Landau theory of a 
second-order phase transition with {*&(z, 9)} as the or- 
der parameter. 



V. CONNECTION WITH STRING THEORY 

In this section, we discuss the similarities and differ- 
ences with superstring theory. For some textbooks on the 
subject, we refer to Refs. [5ll |5g, |57| • In string theory, 
one usually starts with the Polyakov action J53| , describ- 
ing the coordinates X^(a, r), with fx = 0, 1, ..., D — 1, of 
the string propagating in a D-dimensional curved space- 
time with metric G flu (X), 

S = -- f d 2 aVhh a0 d a X' i d0X v G l , v (X) . (91) 

Here a a = {a, r} are coordinates on the worldsheet 
sweeped out by the string, r is the worldsheet time, and 
a runs longitudinally over the string. Furthermore, T is 
the string tension, and h a p is a two-dimensional metric 
on the worldsheet with h = —Dct[h a p]. In agreement 
wit the standard practice in high-energy physics we are 
momentarily using units such that h = c = 1 . We restore 
units when we come to the precise connection with our 
ultracold superstring. In fully quantized string theory, 
one also performs a path integral over these metrics, and 
this leads to the string loop expansion where one sums 
over all two-dimensional surfaces containing an arbitrary 
number of holes. In our setup, the worldsheet of the 
string is completely fixed, and contains no holes, i.e., it 
is just the two-dimensional plane. On the plane, we can 
then make use of the local symmetries of the Polyakov ac- 
tion, that are the reparameterizations of the worldsheet 
coordinates and the Weyl rescalings of the metric. Doing 
so, we can make the gauge choice 



h a p — 



1 
-1 



(92) 



This gauge choice is referred to as the conformal gauge. 
The space-time in which the string propagates is coor- 
dinatized by X^, fj, = 0, ...,D — 1. In quantized super- 
string theory one has that D = 10, but at the classical 
level one can have D = 4 as well. We come back to this 
issue below. It is useful to introduce light-cone variables 



x ± = -L(x°±x D - 1 ) 



(93) 



and X\ i = l,...,D-2. Then X ± and X i describe 
the longitudinal and transversal degrees of freedom of the 
string, respectively. String theory has the special feature 
that there are only transversal physical degrees of free- 
dom. This is because string theory has an additional con- 
straint that can be understood as the equation of motion 
of the worldsheet metric h a p. Defining a = r±a, these 
constraints read in conformal gauge 



d±Xi*d±X''G liV (X) = 



(94) 



and are sometimes called the Virasoro constraints. In 
practice, solving the constraints is difficult, but in the 
so-called light-cone gauge 



X + = 2a'p+T, 



(95) 



where a' = (2ttT)^ 1 and p + is the center-of-mass mo- 
mentum in the X + direction, the longitudinal modes X 
can be eliminated explicitly, at least for certain space- 
time metrics G^ v . The light-cone gauge can always be 
taken as a consequence of the residual gauge symme- 
try after the gauge choice of Eq. (|92|) has been imposed 

The implementation of the constraints in Eq. I|94() 
in the quantum theory leads to the critical dimension, 
namely D — 26 for the bosonic string and D = 10 for 
the superstring. In our condensed-matter setup, these 
constraints are not present. There are physical longitu- 
dinal degrees of freedom, so this makes it different from 
the superstring. However, the longitudinal modes are 
suppressed and at the energy scales we are looking at, 
it suffices to study only the transversal degrees. It is in 
this transversal sector that we connect to string theory. 
To make this connection, we have to specify the space- 
time metric G M „ . A class of backgrounds that has been 
intensely studied in the string literature is that of plane 
wave metrics |5|J,|6(|. The simplest of these backgrounds, 
and also the one relevant for our case, is given by 



As 2 = G fll/ {X)dX tl dX v 



(96) 



-2dX+dX' +H(X i )(dX 



+^2 



dlMX 1 



where H(X l ) is a function of the transverse coordinates 
only. In light-cone gauge, the lagrangian for the string 
propagating in this background now becomes 



D-1 



T- 1 C = 



^ 



dX 1 



fdX' 



V da 



- V{X l ) , (97) 
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where V(X) = -2{a'p+) 2 H(X). To derive this result, 
we simply substitute the background in Eq. (|96|) into 
Eq. (|91(1 . and use the light-cone gauge from Eq. (|95H to 
produce the potential V(X) term in the lagrangian [63|. 
Furthermore, this produces a term proportional to X~ 
that is decoupled from the X 1 . Therefore this term can 
be dropped. In fact X~ is fixed by the Virasoro con- 
straints in Eq. I|94|) , so we only need a lagrangian for the 
transverse degrees of freedom. 

One of the remarkable facts of string theory is that 
its conformal symmetry at the quantum level forces the 
metric to satisfy Einstein's equations in general relativity. 
This is the way in which gravity emerges from string the- 
ory. When there are no other background fields present, 
as in our case, Einstein's equations reduce to a single 
constraint on the function H given by 



D-2 



Atf-E 



if 



{dX 



i\2 



H = 



(98) 



In other words, H has to satisfy the Laplace equation 
in the transverse space. This constraint has to be un- 
derstood on an equal footing as the constraint on the 
space-time dimension. They both follow from a consis- 
tent implementation of the conformal symmetry at the 
quantum level. Since we are not taking into account the 
Virasoro constraints in our system, and hence the con- 
formal symmetry, we therefore also ignore the constraint 
in Eq. (|98() . Doing so, we can work with arbitrary poten- 
tials V(X). When we take D = 4, as we shall below, the 
scalar potential depends on two real fields. 

We now include the fermions, and discuss supersym- 
metry. To make a supcrstring we have to add addi- 
tional terms to the lagrangian in Eq. (|91|l containing the 
fermions in such a way that there is supcrsymmetry. We 
can then impose the conformal or light-cone gauges to ar- 
rive at a supersymmetric generalization of the lagrangian 
in Eq. (|97|) . Alternatively, we can directly study su- 
persymmetric extensions of Eq. I|97|) as two-dimensional 
field theories. The general construction of supersymmet- 
ric two-dimensional field theories with scalar potentials 
V{X) was given in Ref. [62J]. Not all potentials lead to 
lagrangians that can be supersymmetrized. For the case 
of minimal supcrsymmetry with two supercharges, some- 
times denoted by (1,1) SUSY, the potential needs to be 
of the following type 



D-2 



V(X*) = J2((d t W) 2 + GUX)) 



(99) 



where W is a real function, di stands for the deriva- 
tive with respect to X 1 , and the quantities Gi(X) satisfy 
diGj + djGi = together with ^\ G l diW = constant. 
The supersymmetric lagrangian can then be written as 

2T~ 1 £ = d a X l d a X l + i^ l 1 a d a ^-V{X) (100) 



(5), 



-W tj (XW^ - W^{X) v>W, 



with 

The supersymmetry variations are 



(5) 



d^. 



SX i 



eV l 



-i'fdaXU - d'We - G i 1 b e, 



(101) 



(102) 
(103) 



and leave the lagrangian invariant, up to a total deriva- 
tive. Here tp l and e are two-component Majorana spinors, 
and in our model we thus have two Majorana spinors. 
The 7-matrices are related to the Pauli matrices as 
7° = (Jy,^ 1 = icr x and 7 s = a z as before. For more 
details on the spinor conventions, see Ref. |62j . 
Examples of supersymmetric models are given by 



Gi = ae vj X J , 



W(r) =$R + jR 3 



(104) 



where R = ^(X 1 ) 2 + (X 2 ) 2 and a, j3, 7 are arbitrary 
parameters. Plugging this into Eq. (|9^|) leads to J64J 



V(R) = {/3 + 3 7 i? 2 



2 r>2 

a R . 



(105) 



Up to an irrelevant additive constant, the coefficients 
a, /?, 7 can be chosen such that the potential is as in 
our condensed-matter setup. Furthermore, we have that 
W\, = ctCji which leads to mass terms for the fermions, 
and supersymmetry variations of the fermions of the form 



5^ 



ae lj Xjj 5 e 



(106) 



This term rotates the fermions into the bosons, just like 
for the ultracold superstring. If we compute Wtj to de- 
termine the interactions between bosons and fermions, it 
produces complicated interaction terms, 

Wij = Sij (I + 3 7 i?) + XiXj (-7^ + f) , (107) 
as a result of the supersymmetry constraints. 



A. Nonrelativistic limit 

To connect to our condensed-matter setup, we have 
to take the nonrelativistic limit in which only particle 
excitations of the two-dimensional field theory survive, 
and the anti-particle excitations are absent. To illustrate 
this procedure, we start with the bosonic part of the la- 
grangian in Eq. (|100f> . based on two real scalar fields. In 
terms of the complex field 

X = X x + iX 2 , X*=X x -iX 2 , R 2 = \X\ 2 , 

(108) 
the lagrangian reads 



1 



T-'C = ^\d T X\ 2 -\d a X\ 2 -V{\X\) 



(109) 
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where we have reinserted the speed of light c in order 
to take the nonrelativistic limit c — > oo below, and we 
further used that the potential is a function of R only 
since this is the case of interest. 

Using now the mass m*, we decompose the complex 
scalar field in terms of positive and negative frequency 
modes 



X(a,r) 



1 



-im c r 



V2m*T 



b(a, 



„jm c r 



a(a, 



(110) 

and call b the particle field and a the antiparticle field. 
Both b and a are complex. We now substitute Eq. 111(J|) 
into Eq. I|1U9|) and send c — > oo. In this limit, the 
lagrangian becomes first order in time derivatives, and 
particles and antiparticles decouple from each other such 
that we can effectively set a = 0. The remaining terms 
in the nonrelativistic limit are 



C = ihb*d T b 



2m 



-\d a b\'-VQb\), (111) 



where we have reinserted the various factors of h. More- 
over, we have absorbed a mass term proportional to \b\ 2 
into the potential. Remind that we have chosen a po- 
tential of the form given in Eq. (|105fl . so this mass term 
can easily be absorbed into a change of the coefficients 
a or /?7. Notice that this lagrangian precisely coincides 
with the bosonic sector of lagrangian of the ultracold su- 
perstring given in Eq. (|43|l . The fermionic sector can be 
obtained in a similar way. 



cise mathematical connection with string theory in 3 + 1 
dimensions was made. 

The discussions in this article were limited to the case 
of a single string. It is left for future investigation to ex- 
tend the analysis to involve more strings. A complication 
in this case is that for parallel vortex lines, supersymme- 
try is not possible, because of the different way vortices 
and fermions interact with each other. A proposal to 
overcome this problem is to study the interaction of two 
supcrstrings that are both in the center of the condensate, 
but are seperated on the z axis. This would correspond 
to merging and splitting of ultracold supcrstrings. 

The typical fermionic number of particles that is 
needed to obtain supersymmetry is typically around 0.1 
per site. This is rather low, both to control and to ob- 
serve. However, the density is rather high and can be 
estimated to be at least 10 13 cm -3 . Moreover, a dis- 
advantage of a higher fermionic atomic density is that 
this makes also the contribution of the kclvon-fcrmion 
hopping interaction more important. It remains to be 
investigated, whether a change of the varous parameters 
can improve on this situation. 

Apart from the other possibilities mentioned in this 
article, it is also possible to gain experimental insight 
in the system by coupling the vortex motion to reso- 
nant quadrupole modes |42j . This gives the possibility 
to measure the kelvon dispersion directly If the system 
is brought out of equilibrium by populating a high-lying 
kelvon mode, it also opens up the exciting possibility 
to study collapse and revival phenomena between the 
bosonic and fermionic modes. 



VI. CONCLUSION 

In this paper we presented a detailed account of the 
conditions under which the ultracold superstring can be 
created. The requirements for the laser parameters and 
the atomic interactions were given. Moreover we payed 
attention to the experimental signatures of supersymme- 
try. The supersymmetry in the problem was investigated 
by studying the appropriate super algebra. Finally, a pre- 
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